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Abstract: 

Cadabra is a new computer algebra system designed specifically for the solution of problems 
encountered in field theory. It has extensive functionality for tensor polynomial simplification 
taking care of Bianchi and Schouten identities, for fermions and anti-commuting variables, 
Clifford algebras and Fierz transformations, implicit coordinate dependence, multiple index 
types and many other field theory related concepts. The input format is a subset of Tf<]X 
and thus easy to learn. Both a command-line and a graphical interface are available. The 
present paper is an introduction to the program using several concrete problems from gravity, 
supergravity and quantum field theory. 



This paper describes the capabilities of the Cadabra computer algebra system; the software itself, including source code, 
can be obtained from http://www.aei.mpg.de/~peekas/cadabra/ where further documentation can also be found. 
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1 Introduction 

The currently available spectrum of symbolic computer algebra software exhibits, from the per- 
spective of a high-energy physicist, an unfortunate deficiency. On the one hand, problems dealing 
with differential equations, symbolic matrix algebra, special functions, series expansions, poly- 
nomial algebra and so on are adequately covered, often by several different systems. It is not 
too hard to translate a give physics problem in these classes from paper to computer and back. 
Other common problems, on the other hand, such as applying symmetry transformations to La- 
grangians, computing Poisson brackets, deriving field equations, canonicalising tensor expressions 
or performing Fierz transformations, are much harder to handle using existing symbolic computer 
algebra tools. For these classes of problems, which one could label as "field theory problems", it 
typically takes much more effort to convert them from paper and solve them on the computer. 

Granted, there are many smaller systems, often built on top of general purpose systems such 
as Mathematica or Maple, which solve one or more of these problems!]] However, writing such 
packages, or even adapting them for slightly more complicated situations than for which they were 
intended, is often far from easy, and requires a substantial knowledge of the underlying computer 
algebra system. Moreover, combining them, such that e.g. a tensor manipulating package can 
deal with fermions too, is typically a tedious exercise, if possible at all. There are clear technical 
reasons for this, rooted in the design of general purpose computer algebra systems. However, 
rather than dwelling on these deficiencies, it is more productive to actually come up with a 
better solution. 

The present paper is an introduction to a new computer algebra system, called "Cadabra", 
which is designed from the start with field theory problems in mind. By discussing a number 



Among the more often-used recently written packages there is GAMMA for gamma matrix algebra 1], 
grassmann.m for handling of anti-commuting variables 01 and the extensive xPerm for tensor manipulation 3|. 
This list is far from complete and only meant to illustrate how much work one has to do in order to teach a 
general-purpose symbolic computer algebra system about even the simplest concepts which occur in field theory. 
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of explicit sample problems representative of the "field theory" class vaguely defined above, this 
paper is intended as a guide for users rather than computer scientists (for those interested, a 
description of the technical implementation and code design goals has appeared in 0]). 

In which sense does Cadabra differ from other computer algebra systems? The first and most 
easily visible feature is that all expression input is in the form of a subset of Tj]X. Tensor indices, 
Dirac conjugation, derivative operators, commutators, fermion products and so on are all written 
just as in TgX. With a little bit of discipline, one can cut-and-paste expressions straight from 
a paper into a Cadabra notebook. The output, similarly, is typeset as T^X would do it, and 
Cadabra notebook files are in fact at the same time also valid TgX files (the program comes with 
a graphical notebook interface, but can also be used from the command line). 

Secondly, Cadabra contains some of the most powerful tensor expression simplification rou- 
tines available at present. Not only does it use simple symmetry or anti-symmetry of tensors or 
tensor indices in order to simplify expressions^], it also takes into account multi-term relations 
such as the Bianchi identity, as well as dimension-dependent relations such as the Schouten iden- 
tity (section 12. 3p . The program handles commuting as well as anti-commuting tensors, allows 
for multiple index sets (section 12. 5p , and knows about the concept of a dummy index so that no 
special wildcard notation is ever needed when handling tensor indices. 

Thirdly, the program knows about many concepts which are common in field theory. It 
handles anti-commuting and non-commuting objects without special notations for their products 
(section I3.3P , it knows about gamma matrix algebra (section 13 . 1 j) , Fierz identities (section I3.2[) , 
Dirac conjugation, vielbeine, flat and curved, covariant and contravariant indices (section I3.3j) . 
implicit dependence of tensors on coordinates, partial and covariant derivatives (section I2.2[) . 
It has extensive facilities for handling of field theory expressions, e.g. dealing with variational 
derivatives. It features a substitution command which correctly handles anti-commuting objects 
and dummy indices and offers a wide variety of pattern matching situations which occur in field 
theory. 

Fourthly, the program has special support that allows the user to write (and keep) expressions 
in any desired form, e.g. to specify the preferred order of objects. The program does not try to do 
anything smart unless it is explicitly told to do so. However, when desired, it is always possible 
to add any arbitrary simplification step to a list of commands which is executed at every step of 
the calculation. In this way one has precise control over the level of verbosity of the intermediate 
results of a computation. 

Finally, the source code of the program is freely available, with extensive documentation on 
how to extend it with new algorithms and new data types (e.g. twistor variables are planned 
for an upcoming release). The program is completely independent of commercial software and 
relies only on a few other open source libraries and programs. Versions for Linux as well as for 
Mac OS X machines can be downloaded from the web site. 

The following sections discuss these characteristic features of Cadabra, illustrated with many 
explicit calculations and several longer examples (see in particular the Riemann tensor polynomial 
problem in section 12.41 the Kaluza-Klein gravity problem in section 12.51 and the fermion Fierz 
problem in section I3T21) . These examples were chosen to illustrate the typical use of Cadabra, in 
particular the way in which it fills a gap in the existing spectrum of computer algebra software. 
The web site offers a growing collection of longer sample calculations for those readers who are 
interested in applications to more advanced problems. The program is under active development 
and information about updates and new features can also be obtained from there. 



2 I am grateful to Jose Martm-Garcfa for allowing me to use his excellent xPerm code 0] for this purpose. 
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2 Bosonic basics 



2.1 Tensors, indices and symmetries 

Before discussing actual calculations, let us start with a few words concerning notation. This 
discussion can be short because, as mentioned in the introduction, mathematical expressions are 
entered essentially as one would enter them in TgX (with a few restrictions to avoid ambiguities, 
which we will discuss as we go along). In order to manipulate expressions, Cadabra often needs 
to know a bit more about properties of tensors or other symbols. Such properties are entered 
using the ' : : ' symbol. A simple example is the declaration of index sets, useful for automatic 
dummy index relabelling. An example will clarify thisjl 



1 { a, b, c, d }:: Indices . 

2 A_{a b> B_{b c}; 

AabBbc, (1) 

3 ^substitute K /.) ( B_{a b} -> C_{a b c} D_{c> ); 

A a b CbcdDd', (2) 



The automatic index relabelling which has taken place in this substitute command was clearly 
only possible because of the property declaration in the first line. Note how the substitute 
command has also figured out that B_{a b} on the left-hand side is equivalent to B_{b c}, 
without any explicit wildcards or patterns. We will see more of this type of field-theory motivated 
logic throughout the paper. 

Indices can be simple letters, as in the example above, but it is also perfectly possible to put 
accents on them. This can be useful for e.g. notations involving chiral spinors. The following 
example illustrates the use of accents on indices. 



1 A_{\dot{a> \dot{b}}: : Antisymmetric. 

2 A_{\dot{b} \dot{a»; 

*u (3) 

3 Ocanonicalise ! (°/ ) ; 

(-1)4*; W 



Here we also see a second usage of property declarations: the construction in the first line declares 
that the tensor is antisymmetric in its indices. The canonicalise command subsequently 
writes the input in a canonical form, which in this trivial example simply means that the indices 
gets sorted in alphabetical order. Continuing the example above, one can also use subscripts or 
superscripts on indices, as in the example below. 

3 The input and output shown in this paper are essentially identical to those in the graphical interface which 
comes with Cadabra, the only difference being the added line numbers and some minor typographical modifications. 
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1 { a_{l}, a_{2>, a_{3}, a_{4} }:: Indices (vector) . 

2 V_{a_{l» W_{a_{l»: 

3 Osubstitute ! (7„) ( V_{a_{2}} -> M_{a_{2> a_{l» N_{a_{l» ); 

M aia2 N a2 W ai ; (5) 



As is clear from these examples, commands always start with a '@' symbol. The ' ! (%) ' bit 
attached to the commands roughly means that the command should be applied at all levels (the 
' ! ' mark) on the previous expression (the sign). Input lines normally end with ' ; ', but output 
can be suppressed by using ':' instead (as in the last example). 

A guiding principle in Cadabra is that nothing ever has to be declared unless this is absolutely 
needed. This is in contrast to many other systems, where for instance one has to declare manifolds 
and index sets and so on before one can even enter a tensor expression. This makes code hard to 
read, but more importantly, such additional declarations are hard to remember. As an example 
of how Cadabra works, one can e.g. input tensor expressions and perform substitution on them, 
without ever declaring the symbols used for indices. Only when the program needs to generate 
new dummy indices does one need to declare index sets, so that dummy indices can be taken 
out of the right set. The general guideline is that "one only needs to tell the program about the 
meaning of symbols when this is actually required to do the manipulation correctly" . 

In order to declare symmetries of tensors, it is possible to use simple shorthands like the 
Antisymmetric in the example above. More generally, symmetries can be declared using a generic 
Young tableau notation. An object with the symmetries of a Riemann tensor, for example, can 
be declared as in the following example. 



1 R_{a b c d}: :TableauSymmetry(shape={2,2}, indices={0 ,2 , 1 ,3}) . 

2 R_{a bed} R_{d cab}: 

3 Ocanonicalise ! (7.) ; 

( — 1) Rabcd Rabcd] (6) 



The first line indicates that the tensor has the symmetries of the |t |rf | tableau (the numbers in 
the indices argument refer to the index positions). The Ocanonicalise command writes the 
input in canonical form with respect to mono-term symmetries (anti-symmetry in the two index 
pairs and symmetry under pair exchange). The Ricci symmetry is also encoded in the Young 
tableau and will be discussed later. Many tensor symmetries have shorthand notations, so one 
often does not have spell out the tableau (e.g. RiemannTensor or Satisf iesBianchi). 



2.2 Derivatives and dependencies 

There are relatively few surprises when it comes to derivatives. It is possible to write derivatives 
with respect to coordinates, i.e. write d x where x is a coordinate, but it is also possible to 
use indices, as in di with i being a vector index. It is also possible to make objects implicitly 
dependent on a derivative operator, so that one can write dA without an explicit specification of 
the coordinate which is involved here. 
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In order to make this possible, however, derivative objects have to be declared just like indices, 
otherwise the system does not know which symbol (d, D, d, V, . . . ) one wants to use for them. 
Implicit dependencies of objects on coordinates associated to derivatives is indicated through a 
Depends property. Here is an example to illustrate all this: 



1 \nabla{#} :: Derivative . 

2 \partial{#} : : PartialDerivative . 

3 A_{m n} :: Antisymmetric . 

4 V_{m} :: Depends (\nabla) . 

5 

6 \partial_{m p}( A_{q r} V_{n> ) A~{p m}; 

d m p(Aq r V n )A prn ; (7) 



7 Ocanonicalise ! (°/ ) ; 
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0; (8) 

\nabla_{m p}( A_{q r} V_{n> ) A~{p m} ; 
Ocanonicalise ! (°/ ) ; 



(-l)V mp (A gr V n )A m P; (9) 
Ounwrap ! (7„) ; 

(-l)A qr V mp V n A m P; (10) 



Note how the symmetry of a double partial derivative has automatically been taken into account 
(it is part of the PartialDerivative property). This is called "property inheritance". 

2.3 Bianchi, Ricci and Schouten identities 

So far we have seen several examples of so-called "mono-term" canonicalisation, in which sim- 
ple symmetries of tensors are used, which relate one particular term to one particular other 
term. More complicated symmetries are Ricci or Bianchi identities, which relate more than two 
terms ("multi-term" symmetries). Such identities can be taken into account using Young tableau 
projectors. Here is an example to show how this works. 



1 {m,n,p,q,r,s , t#} : : Indices (vector) . 

2 \nabla{#} :: Derivative . 

3 R_{m n p q} : : RiemannTensor . 

4 \nabla_{m}{R_{p q r s}} : : Satisf iesBianchi . 
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The last line is a shorthand, but internally does nothing more than to associate a particular 
Young tableau symmetry to the given tensorQ Here it effectively associates the V operator 
to the Riemann tensor. We can see this in action, for instance, by verifying that the Bianchi 
identity indeed holds, El 

5 \nabla_{m}{R_{p q r s}} + \nabla_{pMR_{q m r s}} + \nabla_{q}{R_{m p r s}}: 

6 @young_pro j ect_tensor ! 2 (°/ ) {ModuloMonoterm} : 

7 @collect_terms ! (%) ; 

(11) 



As expected, by Young-projecting the expression, the Bianchi identity becomes manifest 0]. 

A similar logic can also be used 0] to take into account dimension-dependent identities, 
more generally known as Schouten identities or Lovelock identities. Take as an example two 
anti-symmetric tensors A mnp and B mnp . Then we have 

Amnp-B m nq -^-mnqBmnp — in d — 3, (12) 

which is conventionally proved by using the fact that anti-symmetrising in four indices yields zero 
in three dimensions. Cadabra instead uses Young tableau product rules to accomplish the same 
result 0]. 



1 { m, n, p, q }:: Indices (vector) . 

2 { A_{m n p} , B_{m n p} } : : Antisymmetric . 

3 A_{m n p} B_{m n q} - A_{m n q} B_{m n p}; 

A m npB m nq A mn qB rnn p, (13) 

In four dimensions or higher, this product cannot be simplified any further. Decomposing 
the product using Young projectors and canonicalising with respect to monoterm symmetries 
gives back the input, 

4 { m, n, p, q }:: Integer (1 .. 4) . 

5 @decompose_product ! (%) : 

6 Ocanonicalise ! (%) : 

7 @collect_terms ! (%) ; 

ApmnBqmn Aq mn Bp mn) (14) 

In three dimensions, however, the expression vanishes by virtue of the Schouten identity. 
The same procedure as above now indeed produces a zero, 



4 The alternative is \nabla_{m}{R_{p q r s}} : :TableauSymmetry ( shape={3,2}, indices={l,3,0,2,4} ). 
which makes use of the more general tableau symmetry property which we have seen earlier. 

The "2" following the exclamation mark on @young_project_tensor ! 2(7.) indicates that this command should 
only be applied at "level 2" of the expression, i.e. at the level of the terms of the sum. 
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8 { m, n, p, q }:: Integer (1 .. 3) . 

9 @decompose_product ! (°/ ) : 
10 Ocanonicalise ! (°/ ) : 

n @collect_terms ! (%) ; 

0; (15) 



Some more examples can be found in [4| and in the reference manual. 

As an added bonus, the Young tableau routines which are used internally can also be used 
directly. There are properties both for filled and unfilled tableaux. The following example shows 
an example of the latter, together with the Littlewood-Richardson algorithm and the use of the 
hook formula to compute dimensions of representations, 



1 \tableau{#} : : FilledTableau(dimension=10) 

2 \tableau{0,0Ml , 1} \tableau{a,a}{b,b}: 

3 ®lr_tensor ! (°/ ) ; 
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4 Otabdimension ! (7.) ; 

70785 + 235950 + 154440 + 101640 + 103950 + 13860: 
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2.4 A Riemann tensor polynomial example 

Let us now discuss a somewhat more complicated example of bosonic tensor algebra, which shows 
not only how symmetries of tensors are used, but also how repeated derivative operators work 
and how more complicated canonicalisation problems are solved. 

To this extent we will discuss a certain identity, a proof of which can be found in appendix A 
of 0]. Concretely, the problem concerns two particular third-order polynomials in the Weyl 

tensor Cmnpqi 

given by 

rp (~t mkl (~t (~< pmq , mkl/n /~i pq f~* i~tkmpq/~<l 

ij — * ^jpkq^l ~r a ^jmpq^kl cS^ikjV^ ^ mpq i 

. (18) 

rp /~i (~impqn/~i jk /~ipqmn/~<jk 

& — l ~ / jmnfc < - / ^p q ~r jkmn^ ^ pq ■ 

When evaluated on an Einstein space, these polynomials are supposed to satisfy the identity 

1, 

6 

Proving this is a tedious exercise with Bianchi identities when done by hand. With Cadabra the 
proof is straightforward; the notebook is depicted below. 

The key elements in this example are the declaration of the Weyl tensor as well as the covariant 
derivative, together with the 4th line which links the two, by stating that the covariant derivative 
of the Weyl tensor satisfies the Bianchi identity. 



V l V :j E ij - -V l V i E = 0. (19) 



1 {i, j ,m,n,k,p,q,l,r,r#}: : Indices(vector) . 

2 C_{m n p q}: :WeylTensor. 

3 \nabla{#} :: Derivative . 

4 \nabla_{r}{ C_{m n p q} } : : Satisf iesBianchi . 

6 Eij:=- C_{i m k 1} C_{j p k q} C_{1 p m q} + 1/4 C_{i m k 1} C_{j m p q} C_{k 1 p q} 

7 - 1/2 C_{i k j 1} C_{k m p q} C_{1 m p q} : 

8 

9 E:= C_{j m n k} C_{m p q n} C_{p j k q} + 1/2 C_{j k m n> C_{p q m n} C_{j kpq}: 

10 

n exp:= \nabla_{i}-[\nabla_{j}{ @(Eij) » - 1/6 \nabla_{i>{\nabla_{i>{ 0(E) »: 

We now need to apply (twice) the Leibniz rule to expand the derivatives, and then sort the 
tensors and write the result in canonical form with respect to mono-term symmetries, 

12 Odistribute ! (7.) : Oprodrule ! ('/.) : 

13 Odistribute ! (7.) : Oprodrule ! ('/,) : 

14 

is Qprodsort ! (7o) : Ocanonicalise ! (7o) : @rename_dummies ! (%) : 

16 @collect_terms ! (7.) : 

Because the identity which we intend to prove is only supposed to hold on Einstein spaces, 
we set the divergence of the Weyl tensor to zero, 

17 Osubstitute ! (7.) ( \nabla_{i}{C_{k i 1 m}} -> ) , \nabla_{i}{C_{k m 1 i}} -> ) ; 

CijmnCikmp^ jC n kpq Cymn^ k^ipmq^l pCjq n k 2Cjj' mn ^ iCmkpq ^ pC jknq 
CijmnCikpq^ pC jqnk "^-'ijmn^-'ijkp^ q^mCnqkp ^Cymn^ ? kCijpq\/ 7 pCmnkq 
~^Cij mn ^7 %C jkpq^ k^mripq "V/mnC 'jfcpqV ' j\ 7 fcC 1 'mnpq ~\~ ijmnC ikpq^J 'mV j'Cnfcpg 

1 1 1 (20) 

~^Cij mn ^7 iCmkpq^ nCjkpq ~^Cij mn S7 iCjkpq^ mCnkpq ~\~ ijmnC ikpq^/ ' 7 mC n kpq 
ijranC %kmp^ q^ qCjknp Cij mn ^ kCipmq^ kCjpnq ~~^Cij mn CijkpS/ ? q^ ' qCmnkp 

This expression should vanish upon use of the Bianchi identity. By expanding all tensors 
using their Young projectors, this becomes manifest, 



18 
1!) 



@young_project_product ! (7») : 
Osumf latten! (7.) : 



20 @collect_terms ! (7.) ; 

0; (21) 
This proves the identity (|19p . 



The use of Young projector methods also allows for other calculations which are computationally 
intensive when done by hand. An example is the problem of generating a complete basis of 
monomials of Riemann tensors, discussed in 0]. 
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2.5 Index ranges and subspaces, Kaluza-Klein gravity 



We have so far seen rather simple uses of indices, in which only a single vector space was used. 
Cadabra contains functionality to deal with more complicated situations though. In section [3] 
we will discuss the use of multiple index types as well as implicit indices. Here we will discuss 
another index feature, namely the possibility to "split" indices into two or more subspaces. Let 
us first illustrate this with a simple example. 



1 {M, N, P}: : Indices (space) . 

2 {m, n, p} :: Indices (subspacel) . 

3 {a, b, c} :: Indices (subspace2) . 

4 

s A_{M N> B_{N P}; 

6 @split_index! C/.HM, m, a}; 

A-Mm B m p + Amo, B a p ; (22) 



The first three lines declare three types of indices, labelled as "space", "subspacel" and "sub- 
space2" respectively. The last line of the input performs the split of the dummy index into the 
two subspaces. Instead of using two subspaces labelled by indices, it is also possible to use a one- 
dimensional subspace. An example of how this works can be found in the Kaluza-Klein problem 
discussed below. 

The index split functionality is particularly useful for Kaluza-Klein type problems, in which 
the metric is decomposed according to 

_ ( i^ 1 hmn + 4> A m A n (j) A m 
9lM/ ~{ <t>A n <j> 

(where the usual conventions in four space-time dimensions were used). It is a somewhat tedious 
exercise to compute the Riemann tensor components for this particular metric ansatz. One finds 
for instance that 

RmAnA = ~ V m d n <f> ~ ^ m 05„0 <fT X + ^d^d^' 1 h m nh pq + ^F mp F nq (p 3 h pq , (24) 

where F mn is the field strength of A m . The following notebook computes this expression using 
Cadabra and the index splitting functionality. All intermediate output is suppressed as it tends 
to get rather lengthy. 




1 {\mu,\nu,\rho,\sigma,\kappa,\lambda,\eta,\chi#}: . i n di ce s (full , position=f ixed) . 

2 {m,n,p,q,r,s,t,u,v,m#}: : Indices (subspace, position=f ixed, parent=full) . 

Note the appearance of parent=f ull. This indicates that the indices in the second set span 
a subspace of the indices in the first set. The remaining declarations are standard, 

3 \partial{#}: : PartialDerivative . 

4 g_{\mu\nu} :: Metric . 

5 g~{\mu\nu} : : InverseMetric . 

6 g_{\mu? \nu?} :: Symmetric . 

7 g~{\mu? \nu?} :: Symmetric . 
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8 h_{m n}:: Metric. 

9 h~{m n} : : InverseMetric . 

10 \delta~{\mu?}_{\nu?}- : : KroneckerDelta. 

n \delta_{\mu?}~{\nu?} : : KroneckerDelta. 

12 F_{m n} :: Antisymmetric . 

We will want to expand the Riemann tensor in terms of the metric. The following two 
substitution rules do the conversion from Riemann tensor to Christoffel symbol and from 
Christoffel symbol to metricjj Index patterns like \lambda? match both four- and three- 
dimensional indices. 

13 RtoG:= R~{\lambda?}_{\mu?\nu?\kappa?} -> 

14 - \partial_{\kappa?}{ \Gamma~{\lambda?}_{\mu?\nu?} } 

15 + \partial_{\nu?}{ \Gamma~-C\lambda?}_{\mu?\kappa?} } 

16 - \Gamma~{\eta}_{\mu?\nu?} \Gamma~{\lambda?}_{\kappa?\eta} 

17 + \Gamma~{\eta}_{\mu?\kappa?} \Gamma~{\lambda?}_{\nu?\eta} : 

18 

19 Gtog:= \Gamma~{\lambda?}_{\mu?\nu?} -> 

20 (1/2) * g~{\lambda?\kappa} ( 

21 \partial_{\nu?}{ g_{\kappa\mu?} } + \partial_-C\mu7Mg_-C\kappa\nu.?} } 

22 - \partial_{\kappa}-[ g_{\mu?\nu?} } ) : 

Now input the Rmini component and do the substitution. After each substitution, we 
distribute products over sums. We also apply the product rule to distribute derivatives over 
factors in a product. 

23 todo:= g_{ml m} R~{ml}_{4 n 4} + g_{4 m} R~{4}_{4 n 4} ; 

24 ©substitute ! (7.) ( @(RtoG) ); 

25 ©substitute ! (7.) ( @(Gtog) ); 

26 ©distribute ! (7.) ; 

27 ©prodrule ! (°/ ) ; 

28 ©distribute ! (7 ) ; 

29 ©prodsort ! (7») ; 

We now split the fi index into a m part and the remaining 4 direction (the ! ! version of 
the command makes it apply until the result no longer changes). After that, we remove x 4 
derivatives of the gauge field and write the expression in canonical form, 

30 @split_index! ! C/.){\mu,ml,4}; 

31 Ocanonicalise ! (*/,) : 

32 ©substitute ! (7.) ( \partial_{4}{A??> -> ) ; 

33 Osubstitute ! (7.) ( \partial_{4 m?}{A??> -> ) ; 

In the next step, we insert the metric ansatz (|23p and simplify the result as much as possible. 

34 Osubstitute ! (7.) ( g_{4 4} -> \phi ) : 



35 



Osubstitute ! (7.) ( g_{4 m> -> \phi A_{m> ) : 



36 Osubstitute ! (7.) ( g_{m n} -> \phi**{-l} h_{m n> + \phi A_{m} A_{n> ); 

37 Osubstitute ! ('/,) ( g~{4 4} -> \phi**{-l> + \phi A_{m> h~{m n} A_{n> ): 

38 Osubstitute ! (7.) ( g~{4 m} -> - \phi h"{m n> A_{n> ): 



6 



A library with common substitution rules like these is in preparation. 
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39 ©substitute ! (7) ( g~{m n} -> \phi h~{m n} ): 

40 ©distribute ! (7.) : 

41 ©prodrule ! (°/ ) : 

42 ©distribute ! CD : 

43 ©prodrule ! (°/ ) ; 

44 ©distribute ! CD ; 

45 ©canonicalise ! CD : 

46 ©substitute !! (°/„) ( h_{ml m2> h~{m3 m2> -> \delta_{ml}~{m3} ); 

47 @eliminate_kr ! CD : 



Some derivatives have to be rewritten to a canonical form, 



48 Osubstitute ! (°/„) ( \partial_{mH\plii**{-l}} -> -\phi**{-2> \partial_{m}{\phi> ): 

49 @collect_f actors ! CD : 
bo ©prodsort ! CD : 

si Osubstitute ! CD ( \partial_{p}{h~{n m}} h_{q m} -> - \partial_{p}{h_{q m}} h~{n m} ); 

52 Ocanonicalise ! CD : 

53 ©substitute ! ("/,) ( h_{ml m2> h~{m3 m2> -> \delta_{ml}~{m3} ); 

54 @eliminate_kr ! (7„) ; 

Finally, we replace the derivative of the gauge field with the field strength, 

55 Osubstitute ! ('/,) ( \partial_{n}{A_{m» -> l/2*\partial_{n}{A_{m» + l/2*F_{n m> 

56 + l/2*\partial_{m}{A_{n» ) : 

57 ©distribute ! CD : 
5s ©prodsort ! CD : 

59 ©canonicalise ! CD '■ 

60 @rename_dummies ! CD : 

61 @collect_terms ! CD ', 

- \d m §d n 4>§- 1 + \d p 4>d n h mq h pq - -d mn (f> + ^F mp F nq (f) 3 h m 

1 4 ! 1 ( 25 ) 

+ -d p (j)dq(j)(j) 1 h mn h pq - -d p (t)d q h mn h pq + -dp(j)d m h nq h m - 1 

This is indeed equivalent to (|24p upon writing out the covariant derivative in the first term 
of that equation. 



If required, some of these calculations can be done with fewer lines of input by adding a number 
of default simplification rules; an example of such default rules will be discussed in the next 
section. We will end here the discussion of purely bosonic problems and move on to fermions and 
anti-commuting tensors. 



3 Fermions, Dirac algebra and Fierz transformations 
3.1 Simple gamma matrix algebra 

Cadabra has built-in algorithms for the manipulation of anti-commuting objects, spinors and 
gamma matrices in any dimension. Combined with the option of "suppressing" indices (in our 
examples below we will suppress spinor indices), it becomes possible to write calculations in a 
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natural and compact way. All anti-commuting and fermionic objects as usual need to be declared 
by attaching the appropriate properties to them; we will see many examples of this. 

Let us start, however, with some simple gamma matrix algebra. As an example, we will 
expand the product T sr T r iTkmT ms in arbitrary dimensions in terms of the irreducible T k i and 5 k i 
components. 



1 : :PostDef aultRules( OOprodsort ! (°/„) , @@eliminate_kr ! CD , 

2 OOcanonicalise ! CD , O0collect_terms ! )■ 

This single line adds some default simplification rules to the system, so that we do not have 
to type these by hand at every stage of the calculation. 

We now declare the vector indices, their range, and the symbols used for gamma matrices 
and Kronecker deltas. 

3 {s,r,l,k,m,n}:: Indices (vector) . 

4 {s,r,l,k,m,n}: :lnteger(0. .d-1) . 

5 \Gamma_{#}: : GammaMatrix(metric=\delta) . 
g \delta_{m n} : : KroneckerDelta. 

The declaration for the gamma matrix shows that we are defining an object with implicit 
indices: the spinor indices will be suppressed. The notation _{#} denotes the presence of 
an arbitrary number of indices. 

Next follows the actual calculation. 

7 \Gamma_{s r} \Gamma_{r 1} \Gamma_{k m} \Gamma_{m s}: 
s Ojoin! (°/ ) {expand}: 
9 Ojoin! (7.) {expand}; 



10 ©distribute ! ("/,) ; 

n Ojoin! (7 ) {expand}; 

12 Odistribute ! CD ; 

13 Of actorise ! CD {d} ; 

14 0collect_f actors ! CD ; 



The key ingredient here is the ©join command, which takes two adjacent gamma matrices and 
expands their product in terms of a basis of fully antisymmetrised gamma matrices. In the step 
before (|26p the two commands join the 1st and 2nd gamma matrix, as well as the 3rd and 4th. 
The product is then expanded out and one final gamma join operation is performed. The last 
two lines combine all the ^-dependence in overall factors. 

In the example above, the spinor indices on the gamma matrices were suppressed. The 
GammaMatrix property has turned the gamma symbols into non-commuting objects, which will 
not change order when sorting symbols in a product. It is, however, possible to add the spinor 
indices back in, and use explicit indices. For this purpose, Cadabra knows the concept of an 
"index bracket" , which associates indices to matrix objects like the gamma matrices above, or to 
e.g. vectors or spinors. Here is a somewhat simpler example: 



(-1)((2IV - T ir d + 5 lr d - 5 lr ){2T kr - T kr d + 5 kr - 5 kr d)); 



(26) 



r w (12 - 18d + 8d 2 - d 3 ) + 6 k i(-3 + 6d - 4d 2 + d 3 ); 



(27) 
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1 {a,b , c , d#} :: Indices (spinor) . 

2 \Gamma_{#} : : GammaMatrix(metric=\delta) . 

3 (\Gamma_{m n})_{a b} (\Gamma_{n p})_{b c}; 

4 ©combine ! (7) ; 

(28) 

5 Ojoin! (7.) {expand}: 

6 Ocanonicalise ! (7.) ; 

(J^mp^nn ^mn^np ~i~ ^pn^mn ~\~ $rnp$nn ^mn^np)acj (^9) 



The join algorithm has acted 'inside' the index bracket, on the matrix objects themselves. Con- 
versely, index brackets can be distributed over elements in the sum, and objects can be taken out 
of the index bracket when they are known not to contain implicit indices: 



7 Odistribute ! (7.) : 

8 ©expand ! (7) ; 

(J^mp)ac$nn {^mnjac^np (J^pn^ac^mn {$mp$nn)ac ($mn$np)aci 

(30) 



In this way, component and abstract calculations can be mixed at will. 
3.2 Fierz transformations 

Cadabra can apply Fierz transformations in any dimension to re-order four-fermion terms. As 
an example, consider the following identity for Majorana spinors in eleven dimensions @], 

-e^r^xv^re) = ^^(^ M r ii -" i »V'p)(er r 'r i »"-«re) 

n 

= ^r m ^ p e v] s (sen + io v sm 9e) 



With conventional tools, it would take some time to convert this problem to the computer. 
Proving this with Cadabra, on the other hand, requires little more than defining symbols properly 
and inputting the expression on the left-hand side as one would type it in a paper. Moreover, 
reading off the output is simple as well, because Cadabra produces output which is virtually 
identical to the right-hand side of the equation above. The notebook with comments is displayed 
below. 
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1 {\mu, \nu, \rho} :: Indices (curved, position=f ixed) . 

2 {m,n,p,q,r,s,t#}: : Indices (flat) . 

3 {m,n,p , q,r , s , t#} : : Integer (0 . . 10) . 

4 T~{#{\mu}} :: Antisymmetric . 

5 \psi_{\mu}: : Self Ant i Commuting. 

6 \psi_{\mu} :: Spinor (dimension=ll , type=Majorana) . 

7 \theta: : Spinor (dimension=ll , type=Majorana) . 

8 \epsilon: : Spinor (dimension=ll , type=Majorana) . 

9 {\theta, \epsilon, \psi_{\mu}} : : AntiCommuting 
10 \bar{#} : : DiracBar . 

n \delta~{m n} : : KroneckerDelta. 

12 \Gamma~{#{m}} : : GammaMatrix(metric=\delta) . 

These lines define the properties of all the symbols. We now input (minus) the left-hand 
side of (|3ip . and do a Fierz transformation to bring the 9 and e spinors together, 

13 T~{\mu\nu\rho} e_{\nu}"{s} 

14 \bar{\theta} \Gamma~{r s} \psi_{\rho} 

15 \bar{\psi_{\mu}} \Gamma~{r} \epsilon; 

16 

17 @fierz!(7)( \theta, \epsilon, \psi_{\mu}, \psi_{\nu} ); 

- ^T^f>e s Jr rs T mn T r e^T nm ^ p - -^T^Pe s u 9T rs T mn PT r e^T pn ^ p (32) 

- ^T^e s Jr rs T m ^T r ^ p T qpnm ^ p - ^T^eiev^v^^ e%Ytx qpnm ^ P -, 

This is not yet in the desired form, so we join gamma matrices twice and canonicalise the 
resulting expression, 
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j oin ! (7) {expand} : 

19 Odistribute ! CD : 

20 @eliminate_kr ! (7) : 

21 @j oin ! CD {expand}: 

22 Odistribute ! (7) : 

23 @eliminate_kr ! (7) : 

24 @collect_terms ! (%) : 

25 Ocanonicalise ! (7) : 

2C @rename_dummies ! CD : 

27 @collect_terms ! CD ; 



This is indeed equivalent to the right-hand side of (|3ip . 



(33) 



+ _ T ^Pe p m 9T n e^ u T mn ^ p + — m 9T n ™ r e^T mnpqr i> p ; 



Note in particular once more the simplicity of the input on line 13-17. As advertised, very little 
knowledge of the program is needed in order to be able to read and follow a calculation in a 
Cadabra notebook. 
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3.3 Other assorted topics 

So far we have only seen the substitution command in action on bosonic objects. However, the 
substitution command in Cadabra is aware of anti-commuting objects as well, and will take care 
of signs whenever products in the pattern and the expression are not ordered in the same way. 
It will also refuse to match a pattern if two symbols are declared non-commuting, and do not 
appear in the same order in the pattern and in the expression. 

Anti-commutativity comes in two flavours: self-anticommutativity and mutual anti-commuta- 
tivity. The first is used for objects which carry an index: if tpn is declared self-anticommuting, it 
means that VvVv = — VvVv Below is an example to illustrate these concepts and the substitution 
functionality: 



1 \psi_{\mu} : : Self Ant i Commuting. 

2 { \chi, \psi_{\mu} } : : AntiCommuting. 

:i \chi A~-[\mu\nu} \psi_{\mu} \chi \psi_{\nu} ; 

xA^xipu, (34) 
i ©substitute ! (%) ( \psi_{\mu} \psi_{nu} -> B_{\mu\nu} ) ; 

(-l) X A^B^ X ; (35) 

By declaring x an d to be mutually anti-commuting, Cadabra knows that a sign should 
be picked up when doing the substitution. 

In order to illustrate the meaning of Self AntiCommuting, we undo the substitution com- 
mand above, 

i @pop ("/.); 

X A^x^; (36) 
Now we declare A^ v to be symmetric and canonicalise the expression, 

1 A~{\mu\nu} :: Symmetric . 

2 Ocanonicalise ! (7.) ; 

0; (37) 

Because tp^ has been declared Self AntiCommuting, the program knows that ip^ipv is anti- 
symmetric in its two indices, and has used this to simplify the expression. 



Incidentally, the @pop command provides an undo facility which allows one to go back any 
arbitrary number of steps in the calculation (by simply applying it multiple times) . 

To conclude, let us discuss one more calculation which combines some of the functionality of 
Cadabra discussed above, and also shows how to handle variational problems, especially those 
involving fermionic objects. An example is the variation of a Lagrangian under supersymmetry 
transformations. Consider the following sample calculation, 
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: PostDef aultRules ( @@eliminate_kr ! (%) , OOprodsort ! (%) , @@collect_terms ! (7„) ) 



D{#}: : Derivative. 
\bar{#}: :DiracBar. 
\delta{A??>: : Derivative. 
{m,n,p,q,r,s,t,u,m#}: : Indices (flat) . 
{m,n,p,q,r,s,t , u,m#} : : Integer (0 . . 3) . 

{\mu, \nu, \rho , \sigma, \kappa,\lambda, \alpha,\beta} : : Indices (curved, posit ion=f ixed) . 
{\mu, \nu, \rho , \sigma, \kappa,\lambda, \alpha,\beta} : : Integer (0 . . 3) . 

Declaration of bosonic fields, 

e~{m \mu}: :Vielbein. 
e_{m \mu} : : InverseVielbein. 
g~{\mu\nu} : : InverseMetric . 
g_{\mu\nu}: : Metric. 

Declaration of fermionic fields, 

{ \epsilon, \psi_{\mu} , \psi_{\mu\mi} }: :Spinor(dimension=4, type=Majorana) . 

\Gamma_{#{m}} : : GammaMatrix(metric=\delta) . 

{ \psi_{\mu\nu} , \psi_{\mu} , \epsilon }: : AntiCommuting. 

{ \psi_{\mu} , \psi_-C\mu\nu} }:: Self AntiCommuting. 

{ \epsilon, \psi_{\mu}, \psi_{\mu\nu} } : : SortOrder . 

\Gamma_{#}: : Depends (\bar) . 

\psi_{\mu\nu} : : Antisymmetric . 

Input of the Lagrangian and rewriting such that all indices on gamma matrices are flat, 

L:= -1/2 e e~{n \nu} e~{m \mu} R_{\mu\nu n m} 

- 1/2 e \bar{\psi_\mu} \Gamma~{\mu\nu\rho} D_{\nu}{\psi_{\rho}-}- : 
@rewrite_indices ! (%){ \Gamma~{m n p} }{ e~{n \mu} }; 



In the 1.5th order formalism, the supersymmetry transformation rules are given by 
susy:= { e~{n \mu} -> -\bar{\epsilon} \Gamma~m \psi_\nu e~{m \mu} e~{n \nu}, 



Varying under supersymmetry is now a simple matter of using the ©vary command, giving 
it the rules defined above. This will automatically assume infinitesimal variations. 

Ovary! (L) ( @(susy) ) ; 




ee mn e nv _ _ 




(38) 



e -> e \bar{\epsilon} \Gamma~n \psi_\mu e~{n\mu}, 

\psi_\mu -> D_{\mu}{\epsilon} }: 
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+ \ R pvnm ee mp eY p ^ p e pv e np - \ Y mnp D^D v ^ p ee m >"e nv e? p 

- \ Y mnp %D v D p eee mp e nu e pp - ^ Y mnp ^D^ p eY q ^ a ee qo e mp e nu e pp (39) 
+ \ Y mnp iT p B v i) p €lY q ^ a e qp e m(J e nv e pp + 1 Y mnp 4^D u iP p ee mp eY q iP CT e qu e ncT e pp 
+ \T mnp ^D u ^ p ee mp e nv ZT q ^ (J e qp e pa - 

(The rest of the calculation, in which the result is rewritten using partial integration and a 
Fierz identity to show that it vanishes, can be found in a more extensive notebook on N = 1 
supergravity in four dimensions, available from the web site). 

Various aspects of Cadabra are visible here: the use of multiple index types and conversion 
between them, the user-specified sort order of fields which takes into account commutativity 
properties of tensors, the pattern matching routines which automatically deal with dummy index 
relabelling, the readability of input/output and so on. 

4 Conclusions 

We have discussed the capabilities of the new computer algebra system "Cadabra" , by applying 
it to a number of concrete calculations. There are several aspects of this system which make 
it particularly well-suited to solve field-theory problems. Firstly, it uses (a subset of) TpX not 
only for output, but also for input. Compared to other computer algebra systems, this makes it 
much easier to input complicated expressions without errors, as there is hardly any new syntax 
to be learnt. Secondly, the program has built-in facilities for many concepts which occur in field 
theory, like anti-commuting variables, gamma matrix algebra, implicit dependence on coordinates, 
accents, multiple dummy index sets, canonicalisation of tensor expressions and so on. Although 
there exist other systems which share some of the functionality of Cadabra, the implementation 
in Cadabra was designed from scratch so as to make problem solving resemble as close as possible 
the steps one would follow with pencil and paper. 

The program is entirely built on freely distributable software libraries, i.e. it does not make 
use of any proprietary computer algebra system. The system is available in source code as 
well as in binary form for Linux and Mac OS X computers; the reader is referred to the web 
site |http : / / www . aei . mpg . de/ ~peekas/c adabra/ for download and installation instructions as 
well as the address of the mailing list. 

Given the program's scope and size, and the fact that this is only the first public release, there 
of course remains quite some room for improvements and additions. Planned future enhancements 
include e.g. two-component spinors, automatic canonicalisation using Fierz identities, and twistor 
objects. 
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